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In this work, the anisotropic quantum Rabi model with different conpling strengths of the rotating- 
wave and counter-rotating wave terms is studied by using two kinds of extended coherent states 
(ECS). By the first kind of ECS, we can derive a so-called G-function, by which both the regular 
and exceptional solutions can be given. The exceptional solution are just corresponding to the 
crossing points of two energy levels with different parities, so is doubly degenerate. By the second 
kind of ECS, a general scheme for the eigensolutions is derived analytically in a unified way. The 
zero-order approximation is just the adiabatic approximation, and the first-order approximation is 
actually a generalized rotating-wave approximation. The algebraic formulae for the eigensolutions 
are given explicitly in two approximations. The generalized rotating-wave approximations work well 
in a wide range of two different coupling strengths and the qubit detunings. 

PACS numbers: 42.50.Pq, 03.65.Yz, 71.36.-l-c, 72.15.Qm 


I. INTRODUCTION 

Quantum Rabi model (QRM) describes a two-level sys¬ 
tem coimled to a cavity electromagnetic mode (an oscil¬ 
lator) [l|, which can be used to describe the simplest 
matter-light interactions, and has many applications in 
numerous fields ranging from quantum optics, quantum 
information science to condensed matter physics. In con¬ 
ventional quantum optics 0, the rotating-wave (RW) 
terms are kept and the counter-rotating-wave (CRW) 
terms are neglected. So usually the rotating-wave ap¬ 
proximation (RWA) is employed and analytical closed- 
form exact solutions are available. 

Recently, in the circuit QED systems [H , the cou¬ 
pling between the superconducting qubit and the res¬ 
onator can be strengthened by 10%. In this ultrastrong¬ 
coupling regime, the evidence for the breakdown of the 
RWA has been provided by the transmission spectra Q . 
The remarkable Bloch-Siegert shift associated with the 
counter-rotating terms also demonstrates the failure of 
the RWA [3]. So CRW terms could not be omitted, and 
the full QRM should be considered. Although the nu¬ 
merical solutions to the full QRM are extremely easy to 
obtain, the analytical solutions are however highly non¬ 
trivial. The analytical approximate solutions have been 
obtained at different levels, such as weak energy differ¬ 
ence (A) between two-levels (A/w < 0.5, a; is the cavity 
frequency) and deep strong coupling regime , weak 

and intermediate coupling {g/oj < 0.4) m , and whole 
parameter range |12 h 15I|. The analytic exact solutions 
have been also obtained by many groups . Among 

these exact approaches, the eigenvalues are usually (or 
equivalently) determined by zeros of the derived func¬ 
tions. Two continued fraction techniques are formulated 
on the original Fock space [T^ and optimum extended co¬ 
herent states (ECS) [^, where the built-in truncation is 
unavoidable formally. Recently, Braak presented an ana¬ 
lytical exact solution [Tsil using the Bargmann represen¬ 
tation. A so-called G-function has been derived, which 


actually can be written in terms of Heun functions. Al¬ 
though it is not in a closed form, a built-in truncation is 
not needed formally before the practical calculation, in 
contrast with the continued fraction techniques. Alter¬ 
natively, using the method of extended coherent states, 
this solution was recovered in a simpler and physically 
more intuitive way p^ . 

The anisotropic matter-light interacting systems with 
different RW and CRW coupling strengths have been 
studied for a long time, mostly for the theoretical interest 
previously. The quantum chaos has been studied in the 
anisotropic Dicke model [l^. Recently, the well known 
Goldstone and Higgs modes have been demonstrated in 
optical systems with only a few (artificial) atoms in¬ 
side a cavity, which can be described by a few qubit 
QRM [21 L More recently, the study of an anisotropic 
QRM [^, was motivated by the recent experimental 
progress. This model can be mapped onto the model de¬ 
scribing a two-dimensional electron gas with Rashba (an, 
RW coupling relevant) and Dresselhaus {an, CRW cou¬ 
pling dependent) spin-orbit couplings subject to a per¬ 
pendicular magnetic field [111- These couplings can be 
tuned by an applied electric and magnetic field, allow¬ 
ing the exploration of the whole parameter space of the 
model. This model can directly emerge in both cavity 
QED [2^ and circuit QED Q. For example in Ref. [Ig 
a realization of the anisotropic QRM based on resonant 
Raman transitions in an atom interacting with a high 
finesse optical cavity mode is proposed. 

The exact solutions for the anisotropic QRM have been 
obtained using the Bargmann representatio n [2^ . 
The G-function was obtained by Xie et al. [^, where 
both the regular and exceptional eigenvalues can be ob¬ 
tained. The isolated exact solutions at the level crossing 
was found by Tomka et al. [2^. On the other hand, the 
approximate analytic solutions with explicit expressions 
in a wide parameter regime are not given in literature, 
to the best of our knowledge. 

The paper is organized as follows. In Sec. II, we de- 
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scribe the model of the anisotropic QRM. In Sec. Ill, 
by using ECS technique, we derive a new G-function 
to the anisotropic QRM resembling the compact one in 
the isotropic model, giving not only the exact regular 
spectra, but also exceptional solutions right at all the 
level crossing points. In Sec. IV, by another ECS ap¬ 
proach, we present a generalized rotating-wave approxi¬ 
mation (GRWA) to the anisotropic QRM, the formulae 
for the eigenenergies and the eigenstates are explicitly 
given. A brief summary will be presented finally. 


where A is qubit energy difference, (a) is the photonic 
creation (annihilation) operator of the single-mode cav¬ 
ity with frequency w, 51 and 32 are the RW and CRW 
coupling constants respectively, and (Tk(k = x, y, z) are 
the Pauli matrices. Set r = 32/51 as the anisotropic pa¬ 
rameter. 


III. ANALYTICAL EXACT SOLUTIONS 
WITHIN G-FUNCTION TECHNIQUE 


H. MODEL 


Employing the following transformation 


The Hamiltonian of the anisotropic QRM can be de¬ 
scribed as follows [H, [ 2 ^ 

H = 2 ^0-2 + + 5 i (aV- + aa+) + 32 (aV+ -|- acr_) , 

( 1 ) 
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( 2 ) 

we have the Hamiltonian in the matrix (in units oi h = 

LO = 1) 


Hi = PHP-^ 


( a'^a-b/3 (a-I-a'*')-I- 

V -iA + ^at 



a^a — /3 (a -b a^) — 


( 3 ) 


where A± = (31 ± 3 I) /2 and /3 = ^ 3132 . 

We introduce two displaced bosonic operators with op¬ 
posite displacements 

A\=a^ + /3; Al=a^ - p. (4) 


where D{f]) = exp (/3a’^ — /3a) is the unitary displace¬ 
ment operator, | 0 ) is original vacuum state, |n)^ and 
\n)A_ are called ECS [^ . 


The bosonic number state in terms of the new photonic 
operators and A\_ are 


W)a. 



D{-P) |0) 



Dm io), 


J 


The Hamiltonian in terms of can be written as 


Hi = 


+ [ml — A+ (—^A -b A_) — ^aI_ 

(—^A — A_) -b ^A^ Aa^+ ~ (/3 + m/P) A+ — 2/3A+ -b 2/3^ -b A+ 


( 5 ) 


The wavefunction can be expressed as the following series Projecting a+ (w| onto the Schrodinger equation yields 
expansion using the ECS the recurrence relations for the coefficients 


|A+) = 


f Er=o^e„|n)A+ A 

\EZo'^-fn\n)Aj ' 


( 6 ) 
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fm = 


- Sm-l + (^A — A_) fm + —fm-1 


m — X 
A 


(—— A_) Bm-l + —Bm-2 + {j^L — 1 + 2/3^ + 2A+ — x) fm-1 — {P + /P) fm-2 

2Pm 

I- 


(7) 

( 8 ) 


where x = X-^. + E (E is the energy). Starting from 
/o = 1 , we can obtain all fm recursively, which will be 
very useful later. 

Considering the conserved parity, by the coefficients in 
Eq. (|ni), the wavefunction can also be expressed in the 
ECS of the A_-space as 


|Al_) 


EZo(-^rV^-fn\n)A-\ 

EZo {-irV^-en\n)A- J- 


(9) 


If both wavefunctions lH) and ([9]) are the true eigenfunc¬ 
tion for a non-degenerate eigenstate with eigenvalue E, 
they should be in principle only different by a complex 
constant r' 

OO OO 

V^e„|n)A = r'(-1)” V^/„|n)A; 

n—0 n—0 

OO OO 

'/n^-fn\n)A = r'(-1)” \/^e„|n).4- (10) 

n—0 n—0 

Left multiplying the original vacuum state (0| to both 
side of the above equations, and eliminating the ratio 
constant r' gives 

OO OO OO OO 

y] e„/3” y] /„/3" /„/?", 

n=0 n=0 n=0 n—0 

where we have used 

A!(0|n)^+ = {-irV^.{0\n)^_ = (11) 



x=E+). 


FIG. 1: (Color online) G(x) curves for the anisotropic QRM 
at A = 0.7, <71 = 0.8, <72 = 0.6. The blue (red) curve denotes 
even (odd) parity. 


isolated exact solutions in the isotropic QRM, which are 
just the Juddian solutions with doubly degenerate eigen¬ 
values. The degenerate eigenstates are excluded in prin¬ 
ciple in the solutions based on the proportionality of Eq. 
m used in the present ECS technique. It naturally 
follows that the Juddian solutions are exceptional ones. 
With the G-function (1121) at hand, we can also discuss 
the Juddian solution readily. The G-function is also 
not analytic in x but has simple poles at x = 0 , 1 , 2 ... 
For special values of model parameters < 71 , 32 ) A, there 
are eigenvalues which do not correspond to zeros of Eq. 
m; these are the exceptional solutions. All exceptional 
eigenvalues are given by the positions of the poles x = n 


then we define the following G-function with the help of 
Eq. © 

OO 

G±(x) =y](/„Te„)^", (12) 

n—0 

where -!-(—) in the left-hand side is corresponding to even 
(odd) parity, all coefficients are determined by Eqs. © 
and (©.If 3 i =32 = 3 , the G-function of the isotropic 
QRM [1^ is readily recovered. 

We plot the G-function for A = 0.7 ,31 = 0 . 8,32 = 0.6 
in Fig. [TJ The zeros reproduce all regular spectra, which 
can be confirmed by the numerical exact solutions. The 
energy spectra for A = 0.7, r = 1/2 and 2 are presented 
in Fig. © 

Exceptional solutions: We link the degenerate states to 
the Juddian solutions [ 2 ^. Koc et al. [ 2 ^ have obtained 


E = n-X+. (13) 

The necessary and sufficient condition for the occurrence 
of the eigenvalue is immediately given by 

(^P -e„_i -I- ^ 2 ^ ~ 

which provides a constraint on the model parameters. 
They occur when the pole of G±{x) at x = n is lifted 
because its numerator in Eq. © vanishes. Note that this 
exceptional eigenvalue belongs to the states with both 
even and odd parities, so it is doubly degenerate, and 
should be at the level crossing points without exceptions. 
For n = 0, X = 0 

(A-a_)/„ = o. 
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FIG. 2: (Color online) The spectra and the isolated excep¬ 
tional solutions for the anisotropic QRM at A = 0.7. r = 1/2 
(left) r = 2 (right). 


SO we at most have one exceptional eigenvalue for n = 0 
at 



1 - 1 

1 1 


c/t 


1 

72 



(19) 


the Hamiltonian (ED becomes 


H2 


U'^HU 

f a^a-I-a (a^-I-a) —-y—7(0^ —a) 
\ —Y + l — 0) o'I'a — a [d^ + a) 


( 20 ) 


where a = {gi + ^ 2 ) /2; 7 = (51 - 52 ) /2. 

Two displaced bosonic operators with opposite dis¬ 
placements, different from those in Sec. Ill, are intro¬ 
duced 



(15) 


It follows that the first excited state and the ground 
state only intersects for the CRW coupling weaker than 
RW coupling. The parity in the lowest energy state will 
change in this case, so the first-order quantum phase 
transitions occur at g) , in sharp contrast with the 
isotropic QRM. From Eq. m, the first level crossing in 
Fig. [5] (left) should occur at g7 = 0.9661, consistent 
with the numerical calculations. 

For n = l,x = 1, the condition for the occurrence of 
the exceptional solutions is 


2 (si + 52 ) “1 + 




1 


= 0. (16) 


9i 92 


-|- a, b\_ = — a. 


The number state for the bosonic particles B+ and i?_ 
are 


\n)B+ = 


K)’ 


■D{-a) |0), |n)3_ = 


M 

\fn\ 


■D{a) |0), 

( 21 ) 

then the Hamiltonian can be expressed with the number 
operators of the particles B+ and B- 


( BlB+-a^ 

(^-#+ 7 ( 4 -^+) BlB_-a^ )' 

( 22 ) 

The wavefunction can be expressed as the following series 
expansion using these ECS 


If gi = 52 = 5, it reduces to 

+ (4g2 - 1) = 0, (17) 

which is exactly the same as that in the isotropic QRM. 

For high order exceptional solutions, the condition for 
the occurrence of the exceptional solutions is rather com¬ 
plicated, but can be numerically estimated straightfor¬ 
wardly from Eg. dm) . All of them are at the level cross¬ 
ing points with open circles in Eig. [5] Finally by the 
iteration in Eq. ED , we can formally write the condition 
for the n-th exceptional solution in terms of all coefficient 

fi<n 


Y.r,f, = o, (18) 

i<n 


\B) = 


J2n7 \n)B^ 

±En=’o\/^(-l)"Cn \n)g_ 


(23) 


where -I- (—) stands for even (odd) parity, and Ntr is the 
truncated number of particles H+ and B-. 

Projecting \m)B+ onto the Schrbdinger equation gives 


Ntr 

{m-a^- E) c™ T (-I)’" ^ Rnt,nCn = 0, (24) 

n—0 

where 

A 

Rm,n — Dm,n 7 \Rm,n+l 5 (25) 

where 


where F^ is only model parameters dependent. 


IV. GENERALIZED ROTATING-WAVE 
APPROXIMATION 


D 


mn —B 


+ 


(m| 



\n)B_ 


= (—1)"* exp(—2a^) 


min[m,n] 

E 

k—0 


n!(2a)™+"-2fe 

(m — fc)!(n — k)\k\ 


To facilitate the approximately analytical study, we use 
a unitary transformation 


Note that Dmn = (2a)" ™ exp(—2a^)LJ):j '"(da^) for 
m < n, Dmn = ^ (-2a)™“”exp(-2a2)L)/“”(4a2) for 
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m>n, and Dm,n = 0 if m < 0 or n < 0. Here L{y) is the 
Laguerre polynomial. Eq. (m can be reduced to that 
in the isotropic QRM [3l| if set 7 = 0, r = 1 {i.e. a = g). 
Next, we will perform the approximation step by step. 

Adiabatic approximations: if Ntr = 0, i.e. the zero- 
order approximation, the eigenfunctions for the quantum 
number m are give by (k = 1,2) 

’ (26) 

and the corresponding eigenvalues are 

E^^'> =m-a^T (27) 

Similar to the isotropic QRM [3l| , the zero-order approx¬ 
imation in this technique is just the adiabatic approxima¬ 
tions. It is also the same as the adiabatic approximation 
in Ref. [2^ derived in an alternative way. Note that 


in the adiabatic approximations, the transition between 
states belonging to different manifolds m is neglected. 

Generalized Rotating-Wave Approximations: Beyond 
the adiabatic approximation, the transition between dif¬ 
ferent manifolds should be considered. We will per¬ 
form a further correction by taking into account the 
transition between states belonging to two manifolds m 
and m + 1. The solutions for main quantum number 
m (= 0,1, 2..) can be obtained by selecting two terms in 
Eq. for each m and to-|- 1. Considering =F (—1)"* = 1 
for the implied parity, we have the following determinant 
in a 2-by-2 block starting with m 

( ^ CL -f Rm,m) _ q 

Rm+l^m (1 ^ C : Rjfi-j-2_^m-l-l^ 

where E = m + x. So we readily obtain the eigenenergies 
for each m (= 0,1,2...) as 


— — 0^^ + — {Rm,m ~ 77m-|-l,m-|-l) + (~1) ^ 


Note that for isotropic QRM, it is exactly the same as 
the GRWA result derived in Ref. [ll, [T^ in the isotropic 
QRM. So the present first-order approximation is also 
termed as the GRWA. 

The lowest energy with even parity is 

= l-a^ + liRo.o-Rip) 

([1 - (Ro.o + i?i.i)]' - 4i?o.ii?i,o^9) 

For a given anisotropic parameters r = 1/2 and r = 2, 
by Eqs. (1^ and Eq. (1^ . we calculate the energy levels 
against a = (51 -I- 52) /2 both in the GRWA and the adia¬ 
batic approximation, which are presented in Fig. |31 The 
exact ones from the Sec. Ill are also presented for com¬ 
parison. Obviously, the GRWA results agree well with 
the exact ones qualitatively in the whole coupling regime. 
The crossing properties are all present in the GRWA. The 
results by the adiabatic approximation obviously deviate 
from the exact ones and become worse with increasing A 
and decreasing a. The reason is that the transitions be¬ 
tween states belonging to different manifolds in the true 
physical process are neglected in the adiabatic approxi¬ 
mation, but the dominate transition from the manifolds 
m and m + 1 are taken into account in the GRWA. From 
the energy spectra, the difference between these two ap¬ 
proximations is not very large, but they are essentially 
different. If the wavefunction is involved in some physi¬ 
cal process, the difference should be remarkable. Further 
corrections to the GRWA will only result in qualitatively 
different results, and will not considered here. 


m,m -Sm+l,m+l)] ‘^Rrn,rn-\-lRrn-\-l,rn^ ■; (28) 



FIG. 3: (Color online) The energy levels as a function of 
coupling constant a = {gi + 32) /2 for different anisotropic 
parameters r = 1/2 (left column) and r = 2 (right column) 
and different qubit splitting A = 0.5 (upper panel) and 1.0 
(down panel). GRWA results are denoted by open circles, 
adiabatic ones by solid triangles, exact ones by solid lines. 
Results for even and odd parities are distinguished by black 
and red lines. 


Due to the counter-rotating wave terms, the eigenfunc¬ 
tions and eigenvalues of the anisotropic QRM present 
an open problem because they are not known in any¬ 
thing like a closed form, even the exact solutions reported 
recentlyjl^ and the present new G-functions. No analyt¬ 
ical explicit expressions for the exact eigenvalues for the 














6 


whole coupling range are available in the literature, to 
the best of our knowledge. The analytical explicit ex¬ 
pressions presented in this paper might be practically 
useful. 


V. CONCLUSIONS 

In this work, we first derive a concise G-function, re¬ 
sembling to the compact one in the isotropic model, for 
the anisotropic QRM by using ECS, then obtain quite ac¬ 
curate approximate analytical solutions by another ECS. 
Zeros of the G-function will yield the regular spectra. 
The isolated exact solutions are given by the exceptional 
solutions to this G-function. The condition for their oc¬ 
currence are also derived in the closed form. The crossing 
points of the energy levels satisfy this condition, similar 
to the single-mode QRM. The present analytic solution is 
well defined mathematically, because of no built-in trun¬ 
cations, thereby allowing a conceptually clear, practically 
feasible treatment to energy spectra and many physical 
processes. The explicit expressions for the eigensolutions 
in the GRWA are also obtained analytically by the an¬ 
other ECS. In a wide coupling regime, the GRWA results 
are very close to the exact ones. 

Interestingly, this work adds the anisotropic QRM to 
a list with a compact G-function like 

G±w = f;/Ji±^^adAA')t„b), (30) 

n=0 n X J 


where /„ is determined recursivelv from /g = 1. For the 
isotropic QRM with one-photon 


= —5. 
2 


g”. For the isotropic QRM with two-photons [19|, T^ = 

-7=^(5n,z, Ln{g) is given by Eqs. (48) and (49) there. 

2 V 1-43^ 

For two-mode QRM [33, T^ = — 6n i and Ln{g) 

4vi-s^ ’ 

is given by Eq. (28) there. In the present model, T^ 
is dependent on model parameters and presents for all 
i < n, and Ln — (v'5i<?2) • The qubit-cavity model 
possessing a compact G-function like Eq. o shares 
the common property. The denominator of the parity 
dependent term, i.e. the second term in Eq. dSni), is 
n—X, so the zeros of its numerator will yield the condition 
for the occurrence of the exceptional solutions: isolated 
doubly degenerate eigenstates with eigenenergy x{E) = 
n, which is very helpful to analyze the structure of the 
energy spectra. This list may be expanded by absorbing 
other related models in the future. 
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